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, We present a detailed study of the X-ray and optical spectra of the luminous 

Q ! Seyfert galaxy Mrk 1393, which revealed variable partial obscuration of the active 

^ ! nucleus. The X-ray spectra obtained by XMM-Newton and Swift show moderate 

^ ! absorption with a column density around 3x10^^ cm~^, consistent with a dust- 
reddening interpretation of the steep Balmer decrement seen in recent optical 
spectra. The X-ray flux in the 0.5 to 2 keV band during the XMM-Newton 

'ES \ observation in 2005 and Swift observation in 2006 was a factor 6 brighter than 

I— I, that of the ROSAT All Sky Survey in 1991. In the past 4 years, the broad Ha line 

^ \ brightened by a factor of 4 accompanied by a decrease in the Balmer decrement. 

^ \ A comparison with literature spectra reveals variations in the dust extinction on 

y—^ \ time scales of several years, suggesting that the obscuring material is very close 

^ \ to the active nucleus. These observations indicate that a dust-to-gas ratio as high 

^ J \ as the Galactic value can be present in moderately thick gas in the vicinity of 

O \ the central engine within a few parsecs. We suggest that the obscuring material 

^ \ may be debris disrupted from the dusty torus. 

^ ! Subject headings: galaxies: Seyfert - X-rays: galaxies - quasars: individual (Mrk 

^ ! 1393) - quasars: emission lines 
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Introduction 



The unification scheme for type I and type II nuclei of Seyfert galaxies has been very 
successful in explaining various observational facts: detection of polarized broad emission 



and / or reflection dominated X-ray spectra in Seyfert 2 galaxies (e.g. 



2000 


), presence 


Jaffe et al. 


2004) 



jxion coiumn aensiiies 
Veron-Cettv & Veron 


ge ( 


Mason et al. 


2006|; 



Jaffe et al.ll2004j ). However, we still do not understand the properties of the torus. Clumpy 



distribution is preferred by recent modeling of high resolution infrared im age of nearby 



Seyfert galaxies and the infrared spectral energy distribution (e.g. lHonigll2008l ). but yet little 



is known about the properties of dust clouds, such as their dust-to-gas ratio. Furthermore, 



nosity ( 


Steffen et a 


. 2003 


Bareer et al. 


2005 


2006; 


Eckart et al. 


20061; ^ 


VVane et al. 


2007 


)• 



Partially obscured active galactic nuclei (AGNs) are ideal targets for exploring some of 
these questions. First, since the transmitted light can be directly measured, the estimation 
of th e reddening/ absorpti on column density is much more reliable than in type II quasars 



[e.g. 



Maiolino et al.ll200ll ). This allows a reliable measurement of the dust-to-gas ratio of 
the obscuring material. The transmitted ligh t also allows us to measure properly the nuclear 
properties (e.g., intrinsic nuclear variability. iDong et al.ll2005l ). and puts better constraints 
on the other parameters (such as the size of the broad emission line region). Second, optical 
and X-ray observations also show that at least some of th ese partially obs cured AGNs display 
transitions from type 1.8/1.9 to 1.2/1.5 or vice-versa (e.g. lGoodrichlll989l ). Two different sce- 
narios are proposed to ex plain these transitions: a continuum lurainosity dependent Balme r 



decrement variation (e.g.. IWamsteker et al.lllQQOl : iTran et al.lll992l : IShapovalova et al.ll2004l ). 



and the dust absorbing material moving in/out of line sight. Clearly, simultaneous X-ray and 
optical observation can give important constraints of these models. If the latter scenario is 
confirmed, these transitions carry important information of the distribution and kinematics 
of absorbing material. 

We initiated a program to obse rve partially obscu red AGNs and quasars in X-rays 
with the XMM-Newton Observatory (jjansen et al.ll200ll ). The targets were selected from 
the Sloan Digital Sky Survey (SDSS, York et al. 2000) by their large Balmer decrements 
( IDong et al.ll2005l ). In this paper, we report of the detections of X-ray absorption and large 
variations of the Balmer decrement in one of these objects, Mrk 1393 {z = 0.0544). The X-ray 
observations were performed by XMM-Newton and Swift. Mrk 1393 appears to be a partially 
obscured Seyfert galaxy due to its Balmer decrement as large as H«/H/3 = 7.7 measured 



from its SDSS spectrum. Mkn 1393 is a Seyfert 1 galaxy with prominent broad H/3 (IKinman 



3 



19831 ). and is also included in the broad line AGN sample of Hao et al. (2005). The reddening 
corrected Ha luminosity of 10^^ erg s~^ makes this AGN a borderline object between Seyferts 
and QSOs. It is also infrared luminous with = —25.03 and was detected by ROSAT in 
X-ray (Anderson et al 2003). With a radio flux of 23 mJy detected in the FIRST survey 
(Becker, White & Helfand 1995), it is formally radio loud with an Rj. = log{fiAGHz/ fr) = 1-42 
(Ivezic et al. 2002). The paper is organized as follows: the X-ray data obtained by XMM- 
Newton and Swift are analyzed in §2. The optical spectroscopic observations and analysis 
are presented in §3. The properties of the obscuring material is discussed in §4. Throughout 
this paper, luminosities are calculated assuming a ACDM cosmology with f2M=0.27, f2A=0.73 
and a Hubble constant of Hq=72 km s~^ Mpc~^, corresponding to a luminosity distance of 
D=236.2 Mpc to the galaxy. 



2. Results from X-ray data 
2.1. XMM-Newton and Swift observations and data reduction 



Mrk 1393 was ob served by XMM-Ne wton on 2005 July 20 with the EPIC PN fIStriider et al 



20011 ) and MOSl+2 flTurner e all 120011 ) for 4423s, 6565s, and 6570s, respectively. The PN 
and the MOS detectors were operated in full-frame mode with the thin filter. Due to some 
strong background flares at the beginning of the observation, part of the data had to be 
excluded from further analysis. After excluding the time intervals with a hard X-ray back- 
ground {E > 10 keV) of more than 10 counts s the net observing times are 4099s, 5324s, 
and 5335s for the PN, MOS-1, and MOS-2, respectively. The XMM-Newton data were 
reduced with XMMSAS 7.1.0. Only single and double events were used for the PN data 
analysis, and single to quadruple events for the MOS data analysis. Data were also selected 
with the quality flag set to FLAG=0. The source counts are extracted from a circular region 
with a radius of 90'.' around the source position. Background counts were accumulated in a 
circular region of the same radius close by. The net source count rates in the 0.2-10 keV band 
are 2.087 cts s~^ for the PA^ detector and 1.577 cts s~^ for M0Sl-h2 detectors, indicating 
that the data are not significantly affected by pileup. The light curves were binned in 100s 
bins. The X-ray spectra were rebinned using grppha version 3.0.0 with at least 20 counts 
per bin. The redistribution and auxiliary response matrices were created by the XMMSAS 
tasks rmfgen and arfgen, respectively. 

During the XM M-Newton observatio n, Mrk 1393 was observed with the on-board Op- 



tical Monitor (OM; iMason et al. I [20011) with the UVM2 filter. The OM data were re- 
duced with the XMMSAS task omichain. The measured UVM2 magnitude from the source 
list is M2=16.17±0.08 mag. Correcting for the Galactic reddening of £'b-v=0.059 mag 
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( ISchlegel. Finkbeiner. fc Davislll998l ) results in M2=15.60 mag. 



Mrk 1393 was observed by Swift (IGehrels et al. Il2004f) on 2005 September 10 for a total 
exposure time of 9.6 ks. The 5*^/1^- XRT ( Burrows et al.l2005 ) observation was performed in 
photon counting mode (IHill et al. 1120041 ). The event file of the observation was created by 
using the Swift analysis tool xrtpipeline version 0.11.4. Source photons were extracted in a 
circle with a radius of 47'.' centered on the source, and the background was selected from 
a source-free region with a radius of 188'.'. Only single to quadruple events in the energy 
range of 0.3-10.0 keV were selected for further analysis. Source and background spectra 
were extracted from the event file by using XSELECT version 2.3. Spectra were rebinned 
within grppha 3.0.0 to have at least 20 photons per bin. The auxiliary response files were 
created by the Swift tool xrtmkarf. We used the response matrix version 010 with a grade 
selection 0-12. The data are not significantly affected by pileup due to the low count rate 
(~0.2 cts s~^). The XRT light curve was created by a program in having at least 

200 photons per bin (see e.g. Nousek et al. 2006 for details). 

In additio n to the X-ray data , we also obtained photometry with the UV/Optical Tele- 
scope (UVOT; iPoole et al. II2008I ) in the UVM2 filter. Source photons were extracted from a 
circular region with r=5", and the background in an annulus around the source with an inner 
radius of 7" and an outer radius of 20'q The UVOT tool uvotsource was used to determine 
the magnitudes and fluxes. We searched for UV variability in the individual segments of the 
observation and did not detect any. The average observed magnitude based on the co-added 
images is muvM2 = 15.88 ± 0.01 mag (corr esponding to a flux of / = 2.11 x 1 0~^^ ergs s~^ 
cm^ corrected for Galactic reddening (jSchlegel. Finkbeiner. fc Davislll998l ). This value 
is consistent with the UVM2 OM magnitude of M2adj = 15.90 ma^, indicating no significant 
change in the UV flux between the XMM-Newton and Swift observations. 



2.2. X-ray Spectral Analysis 



The X-ray spectra are fitted by using XSPEC (lArnaud I Il996l . v.l2.3.1x). Errors of 
derived parameters are quoted at Ay^=2.7. The Gal actic HI column density in the direction 
of the source is 4.8x10^° cm~^ (IKalberla et al.ll2005l ) and is always taken into account in all 



1 MIDAS = ESO's Munich Image Data Analysis System 

^At the SDSS u-band, the radius containing 50% of Petrosian flux is about 2 .0. 



GruDe et al. 


2007; 


GruDe et al. 


2008) 
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the spectral fits as listed in Table 1. 



The XMM-Newton spectra in the 1-10 keV energy band can be fit by a power law model 
with additional intrinsic neutral absorption. The intrinsic absorption column density at the 
resdshift of the AGN is (2.4±0.6) x 10^^ cm~^. The X-ray spectrum in the range 1-10 keV 
is somewhat flat {T = 1.50 ± 0.06), even for typical Seyfert 1 galaxies (e.g. iLeighly Ill999 



Grupe et al. Il200ll ). There is no apparent iron Ka line at the corresponding energy. The 
upper limit to a narrow {a = 0.1 keV) iron Ka line is 30 eV. 

A flat spectrum may r esult either from partial covering absorptiorB or from strong 



reflection (reprocessing) (e.g. iFabian et al.lll989 



Gallo et al.ll2006l : lGrupe et al. 1120081 ). Here 



we examine both possibilities. A reflection model [pexrav in XSPEC) does not improve the flt, 
yielding a best fltted normalization R of the reflection component of zero and an upper limit 
of i? < 1.52 (flxing the inclination to 45° and the metal abundance to the solar value). This is 
consistent with the upper limit derived from the equivalent width of the neutral iron line. An 
excess absorption and a flat photon index are still required, with the same values as the simple 
power law flt (F = 1.44too5 and A'^^ = 2.2to5 10^^ cm~^). Neither does a partial covering 
model improve the flt. The covering fraction (0.4-1.0), though consistent with a full coverage, 
is poorly constrained, and so is the absorption column density (A'^^ = 3.8l2'i x lO^^cm"^). 
The photon index remains flat (F = I.48I007). thus conclude that the excess absorption 
and a flat continuum index inferred from the 1-10 keV spectra are caused by neither partial 
covering nor reflection. 

We now consider the whole 0.1-10 keV energy band. Extrapolating the 1-10 keV best-flt 
power law model down to lower energies under-predicts the soft X-ray flux by a factor of more 
than three below 0.6 keV, suggesting the presence of a possible strong soft X-ray excess. It is 
known that in some cases ionized absorption models can mimic soft excess (e.g., Komossa & 
Meerschweinchen 2000, Done et al. 2007). We flrst consider an ionized absorption model to 
examine this possibility. We flt the simple model ahsori available in XSPEC. This produces 
a marginally acceptable flt with a x^=703 for 652 d.o.f. (corresponding to a probability 
level of Pr — 8%). The absorber has moderate ionization (with an ionization parameter 
^ = L^/AnR'^n ~ I8I3 ^^S cm ) and is considerably thicker {N^ = 1-0 ± 0.1 10^^ cm~^) 
than in the cold absorption model. The flt also yields signiflcantly steep X-ray spectrum 
(F = 1.79 ± 0.03). However, systematic deviations are still seen in the residuals both in 0.5 
to 1.5 keV range and above 5 keV (see[T]). The systematic deviations may be attributed 
partly because the ahsori model does not fully represent the spectrum, and partly to the 



''Partial covering may be caused by a collection of absorbing blobs which are individually smaller than 
the X-ray emitting region, such as broad line emission clouds. 
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multi-phase of the ionized absorber. In order to examine the latter possibility, we add a 
second ahsori to the model. This results in a Ax^ = 46 for two more free parameters, and 
all systematical residuals below 2 keV disappear now. The two components have ionization 
parameters ^low = ^-^^/i and ^high = 2591^^4, ^"^^ column densities bA'^]^^ 10^^ cm~^ and 
2.6^f 10^^ 'low' and 'high' ionization phases, respectively. 

Next, we consider the soft X-ray excess to be a real emission feature rather than induced 
by ionized absorption, and add a soft X-ray component in the spectral fit. Adding soft excess 
emission often causes degeneracy between the strength of soft excess and the amount of 
absorption. Since X-ray absorption is our main concern in this work, we focus the following 
analysis on detecting and measuring the amount of any possible X-ray absorption, by testing 
various spectral models of the soft excess, namely, black-body, ionized plasma, and power- 
law. 

One of the common spectral shapes reproducing the soft X-ray excess emission in AGN is 
a black-body. We first test a power-law plus a black-body model, both absorbed by the same 
amount of (cold) absorption. However, the fit is hardly acceptable {x^/d.o.f = 720/657), 
with a probability level of only ^ 3%. The fitted absorption column density and photon 
index are l.Olg 3 x 10^^ cm~^ and T = IAS'^q'q'^, respectively. We also use a broken power-law 
to account for the spectral steepening at the low energies, and find it unacceptable either, 
with Pr ~ 5%. 

Soft X-ray lines are commonly seen in Seyfert 2 and intermediate Seyfert galaxies, 
which is interpreted as arising from extended photoionized gas, such as the narrow line 
region (NLR). This component is less prominent in Seyfert 1 galaxies because of strong 
nuclear X-ray emission. In the case of Mrk 1393, however, this component may become 
prominent given the likely strong attenuation of the nuclear continuum below 1 keV (see 
above). Here we use a thermal plasma model {meka in XSPEC), which is a reasonable 
approximation of the spectrum of photoionized gas given the spectral resolution and signal- 
to-noise ratio of the current X-ray observations. We leave abundances and temperature as 
free parameters. Acceptable fits can be achieved only when the plasma emission is free 
from the excess absorption required for the power-law continuum. This is consistent with 
the plausible nature of as being extended emission. This model yields an acceptable fit 
with x^/d.o.f = 666/650 {Pr = 32%), and a similar absorption column density Nh = 
S.Otoi X 10^1 cm-2, and a flat continuum T = l.^Stof^, kT = 0.166 ± 0.004 keV, and a 
low abundance of Z/Zq = OmStoZt 

(see Table 1). The low abundance suggests that this 
component is dominated by continuum rather than emission lines. Thus, it raises a big 
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concern for interpreting UT as from an photoionized gas §. The total luminosity of this 
component is ~ 4 x 10"^^ erg s~^ in 0.5-2.0 keV, higher than those of most Seyferts but 
comparable to the luminous extended X-ray emission in NGC 6240 (Komossa et al. 2003). 
This large soft X-ray luminosity we observe is unexpected because the [OIII] to soft X-ray 
flux ratio in Mrk 1 393 is 20 times lower than the typical values found in Seyfert 2 galaxies 



(jCappi et al.lll996l ). So it is likely that this soft X-ray emission component in Mrk 1393 is 
diffuse, possibly contributed by photoionized gas (e.g. the inner NLR) and/or thermal gas, 
or scattered light of the continuum. 

Though an absorbed power-law plus an unabsorbed meka component gives an acceptable 
flt, there are systematic residuals in the soft energy band, however. In fact, we flnd that 
adding a second soft X-ray component with the same absorption column density signiflcantly 
improves the flt at low energies. Using a black body model as the (nuclear) soft X-ray 
component reduces the by Ax^ = 23 for 1 more free parameter (643/649). The parameters 
of the best-fit are listed in Table 1, where the abundance of the meka model is fixed at the solar 
value. The solar abundance also ensures that this component can be interpreted as emission 
from photoionized gas. There are little changes in the column density (2.41q;3 x 10^^ cm~^) 
and the photon index (1.481q q4). the plasma emission (meka) the temperature remains 
unchanged (0.167±0.005 keV), while its luminosity is reduced by ~ 40% compared to the 
previous value (2.6x10^^ erg s~^). Using a broken power-law model instead of the power-law 
plus blackbody model yields a similarly good fit and similar model parameters (see Table 1), 
with a luminosity of 2.2 x 10^^ erg s~^ in 0.5-2.0 keV band for the meka component. 

Based on the above results, extensive fits with various plausible spectral models, we 
find that the amount of 'cold' X-ray absorption is surprisingly robust, independent of the 
choice of spectral models in the soft X-ray band. We conclude that, for Mrk 1393, the excess 
'cold' X-ray absorption is inevitably required, and its measured column density is relatively 
insensitive to the choice of the spectral shape of the soft X-ray excess components. Both the 
soft excess and excess absorption can also be accounted by two-phased ionized absorption 
medium. 

The SWIFT XRT spectrum has a lower signal-to-noise ratio compared to the XMM- 
Newton spectrum. We simply apply those best-fit models from the XMM-Newton data to 
the Swift spectrum. We find that all these models fit the Swift spectrum quite well, yielding 
almost the same model parameters as those derived from the XMM-Newton data. Therefore, 



^Note that a similarly low abundance was obtained in a large numbe r of other sources ; for instance, one 
example is the Seyfert 2 galaxy NGC 5252 with thermal emission model (jCappi et al.lll996f ). This commonly 
seen effect is generally traced back to multiple components or non-equilibrium effects in the emitting plasma. 
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the two spectra are well consistent with each other. In Table 1, we list results of only one of 
the best model fits for demonstration. 

To summarize, the XMM-Newton spectrum requires intrinsic absorption of 'cold' gas 
with a column density A'^^ ^ 2.2 — 3.3 x 10^^ cm~^. The power-law spectrum above 1 keV 
has a photon index of F ~ 1.4 — 1.5, marginally flatter than typical Seyfert 1 galaxies. 
An additional emission component, probably a combination of scattered light and emission 
lines from collisionally or photoionized gas, is inferred to be present with a luminosity of 
a few times 10^^ erg s~^ in 0.5-2.0 keV and an equivalent temperature of 0.16 keV. These 
parameters do not vary much for the models considered here. Alternatively, a two phased 
ionized absorber may account for both the excessive absorption and the soft X-ray excess. 
The amount of 'low' ionization phase medium is a factor of two larger than that in the cold 
absorption model. 



2.3. X-ray variability 

The XMM-Newton observation lasted less than 7 ks, during which no significant X-ray 
variability is found. To search for variability in the Swift observation covering 19 hr in total, 
the count rates are corrected according to the exposure maps to compensate losses induced 
when the source was positioned on bad CCD columns. The light curve (Fig [3]) shows some 
variability on the order of 20% within an hour. The curve is rebined that each bin contains 
100 counts. A x^-test gives = 41.8 for 11 degree of freedoms, which has a null probability 
of less 2xl0~^. Furthermore, we also compare the average Swift fluxes in the 0.3-1.0 and 
2.0-10.0 keV bands with the XMM-Newton measurements. The mean fluxes do not show 
any significant variability on month timescales. 



Mrk 1393 was detected during the ROSAT All Sky Survey (RASS; IVoges et al. 1119991 ) 
with a count rate of 0.033±0.012 cts s~^ in the 0.1-2.4 keV band and a hardness ratio 
of HR1=(H-S)/(H+S)=0.73±0.43 with the H and S being the coun ts in the hard (0.5-2.0 



keV) and the soft (0.1-0.4 keV) bands, respectively (jVoges et al.ll2000l ). The observation was 
made between August 6 and August 8, 1991. For comparison, we estimate the ROSAT count 
rate for the best-fit XMM-Newton model using the ROSAT PSPCC response. The broken 
power-law plus meka model predicts an RASS count-rate of 0.189 cts s~^ in the 0.5-2.0 
keV band. The power-law plus black body soft excess and meka model gives a very similar 
count-rate. With the observed hardness ratio and count-rate in the whole band, we estimate 
the observed RASS count rate in the 0.5-2.0 keV band to be 0.028 cts s~^. Therefore, the 
soft X-ray flux increased by a factor of 6 since the ROSAT observation if we ignore th e 



inter-mission calibration uncertainty, which is estimated to be within 10% fISnowderu |2002| ) . 
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In passing, we note that the variation of intrinsic X-ray luminosity might be larger than 
this apparent value. Thermal plasma emission with similar temperatures has been observed 
in Seyfert 2 galaxies, which is emitted on relatively large scales, and is likely not variable 
on time scale less than ten years (e.g. iBianchi. Guainazzi. fc Chiabergd l2006l ). The meka 
component alone predicts a ROSAT count rate in the band (0.5-2.0) keV of 0.027 cts s~^ 
with 15% uncertainty due to its normalization, very close to the observed RASS count-rate. 
If the meka component comes from the emission from the extended photoionized gas, it 
might be not variable. Thus the actual variability amplitude can be even larger. 



3. Optical Observations and Emission Line Measurements 

3.1. Spectra Obtained at Xinglong Station 

The optical spectroscopic observations were made on September 4-5, 2005, and January 
03-08, 2006 at Xinglong station, Beijing Astronomical Observatory (BAG). Low resolution 
spectra were taken with the OMR spectrograph mounted on the BAO 2.16 m telescope, 
using a Tektronix 1024 x 1024 CCD as a detector. A 300 line mm~^ grating was used, blazed 
at 5200 A for the 2005 September 04-05 and at 6000 A for the 2006 January 03-08 runs, 
respectively. Two exposures of 60 minutes each were taken for each of the six nights at 2005 
September 04-05, 2006 January03, 05, 07, and 08, and one exposure of 60 minutes each for 
each of the two nights at 2006 January 04 and 06. The seeing ranged from ~ 2".0 — 4". 5 
during the observations, and we used slit widths of 2".0 for all of the observations when 
the seeing was less than 3".0, except for the three nights at 2006 January 03-05 when the 
seeing was the largest (~ 4". 5). The resolution of the BAO spectra is ~ 10 A FWHM as 
measured from the night-sky lines. The CCD reductions, including bias subtraction, flat- 
field correction, and cosmic-ray removal, were accomplished with standard procedures using 
IRAF. Wavelength calibration was carried out using Fe-He-Ar lamps taken at the beginning 
and the end of the observations. The accuracy of wavelength calibration was better than 
1 A. Flux calibration was derived with observation of a KPNO standard star, which was 
observed at the beginning and end of each run. 

Spectra obtained in the four nights with the best weather condition (seeing ~ 2".0; 2005 
September 04-05 and 2006 January 07-08) agree with each other within 10%, indicating that 
1) the uncertainty of our flux calibration is better than 10%, and 2) the source is not 
variable on short (intranight) and medium (a few months) time-scales. All of the spectra 
were corrected for the Galactic extinction and combined to increase S/N ratio. The result 
is displayed in FiglH 
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3.2. SDSS and Other Spectra 

The SDSS spectrum of MRK 1393 was taken on 2001 March 22 with an exposure 



time of 3041 s and a spectral resolution of about A/ A A = 1800 flAdelman-McCarthy et al 



20071 ). The seeing during the exposure was around 2'/2. We started with the 1-d spectrum, 
processed by the SDSS pipehne, and corrected it for the Galactic reddening. The spectrum 
is also shown in Fig HI Broad Ha appears much weaker than in the BAO spectrum, and 
the broad H/3, which can be easily identified in BAO spectrum, can be hardly spotted in 
the original SDSS spectrum. The spectrum show prominent stellar absorption lines, such as 
Call, HON, Mgib, NalD, suggesting that the continuum is dominated by starlight. 



A KPNO spectrum taken on 1993 April 24 was kindly provided by Owen f lOwen. Ledlow. fc Keel 



I995I ). The spectral resolution is 5.6 A (FWHM). Owen et al. estimated that the effective 
aperture size of the 1-d spectrum is around 3" in diameter with an uncertainty of less than 
30%. The aperture includes about ~95% of the AGN light (assuming a point-like source). 
The contribution of starlight should be comparable to that in the SDSS spectrum. While 
both spectral resolution and the signal-to-noise ratio of the KPNO spectrum are lower than 
the SDSS spectrum, major stellar absorption lines, such as Call, MgIb, Nal can be easily 
identified (FiglH). 



Two more optical spectra can be found in the literature (jKinmarull983l : iMorris fc Ward 



19881 . hereafter, MW88). The spectrum taken in 1984 by MW88 shows prominent broad 
Balmer lines up to He and a moderately blue continuum, which is very different from the 
SDSS (2001) and the BAO (2005-2006) spectra. There is no observation time in the Kinman's 
paper, though it must be taken before August 1981. These authors also listed the emission 
lines they measured. 



3.3. Variations in the Continuum and Emission Lines 

In order to compare the strengths of emission lines among different observations, we 
re-calibrated all the spectra using the [OIII]A5007A line flux, assuming it is not variable 
and the [OIII] emitting region is not spatially resolved. Constancy of [OIII] line flux over 
a time scale of ten years probably is a good approximation, however, the point-like source 
approximation for NLR is a concern. We noticed that the size (containing 98% line flux) of 
[OIII] emission line region for Mrk 1393 is 1 kpc or 0."9, according to the scaling relation of 
Bennert et al. (2002). This size is smaller than the spatial resolution, but is still comparable 
to the seeing disk and the aperture size used for the spectrum extraction, thus we warn that 
the absolute flux of the broad emission line and continuum flux may have uncertainties of 
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order of ten percent. The [OIII]A5007 emission line flux measured from the stacked BAO 
spectrum is only 83.5% of that measured from the SDSS spectrum. The difference can be 
ascribed most likely to the aperture effect. The average seeing disk during our observations 
was larger than that of the SDSS observation, while the apertures we adopted to extract 
1-d BAO spectra are less than that for the SDSS spectrum (about 2 ".0 x 2 ".5 for BAO vs. 
3" in diameter for SDSS). Assuming [OIII] is not variable over timescales of a few years, 
we multiply the BAO spectrum by a factor of 1.198 to match the [OIII] flux measured in 
the SDSS to spectrum. We note that the [01], [SII] and [Nil] fluxes are consistent with 
such a correction factor. The [OIII] flux measured in the KPNO spectrum is a factor of 3.3 
higher than that in the SDSS spectrum, while the correction appears somewhat smaller in 
the blue part as suggested by [OII]A3727 flux (3.0). This can be explained by either effects 
of differential atmospheric refraction or more extended [Oil] emission region. We divide the 
spectrum by 3.3. 

After recalibrating the fluxes, we convolved the SDSS and KPNO spectra with Gaussians 
to match the instrumental resolution of the BAO spectrum. These spectra are displayed in 
Fig. m It is evident that the strength of both the Ha and H/3 broad lines varied significantly 
between different epochs, and the broad H/3 line varied by an even larger factor, whereas 
all of the narrow emission lines remain constant (note that we used only [OIII]A5007 as 
the calibrator). The continuum shapes and intensities of the three spectra also appear to be 
quite different. The continuum in the KPNO spectrum is only ~ 12% higher than that in the 
SDSS spectrum, while the overall shape of the former is significantly bluer than the latter; 
this is very likely due to different nuclear contributions and/or different nuclear continuum 
slopes. 

The continuum mod eling and emission line measurements were carried out in the way 



as described in detail in IZhou et al.l (120061 ) for the BAO, SDSS and KPNO spectra. The 



host galaxy sta rlight is subtracted using the templates constructed using the ICA method 



(ILu et al.ll2006l ). The nuclear component is described by a reddened power-law. The emission 
lines are modeled with Gaussians: one Gaussian for each of the narrow lines, and two 
Gaussians for each of the broad lines. We assumed that the Ha and H/3 broad lines have the 
same redshifts and profiles. Bad pixels are masked in modeling the continuum and emission 
lines. This yields satisfactory fits to all the data. The SDSS spectrum can be fairly well 
modeled by starlight alone suggesting that the nuclear continuum is weak. In contrast, a 
reddened power-law component is required in the BAO and KPNO spectra. 

The Balmer decrement of the broad line component displayed large variations on time- 
scales of several years. We found, for the broad line component, Ha/H/3 = 5.95ib0.23 for the 
BAO spectrum (taken at the end of 2005 and early 2006), 7.69±0.73 for the SDSS spectrum 
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(2001), and Ha/H/? = 3.33 ± 0.12 for the KPNO spectrum (1993). We also estimated the 
broad hne Balmer decrement to be 3.04 in 1984 using the integrated Ha and H/3 hne fluxes 
given by MW88, assuming a constant narrow Ha, H/5, [NH] hne flux over ten years. This is 
rectifled because aU other narrow hne fluxes quoted in MW88 are in good agreement with our 
measurement in the SDSS spectrum. The flat Balmer decrement is consistent with its blue 
continuum spectrum. The overall variation of the broad Ha flux and the intrinsic reddening 
estimated from the Balmer decrement are plotted in Fig. [51 The reddening i s estimated 



from broad line Balmer decrement, assum ing an intrinsic value of Ha/H/3=3.0 (IDong et al. 



2005; 


Zhou et al. 


2006; 


Done et al. 


2008) 



4. Discussion 



4.1. A Reddening Origin for the Variable Balmer Decrement 

Dong et al. (2008) showed that the Balmer decrements for blue quasars are in a very 
narrow range around 3.0 and do not correlate with either the continuum slope, Eddington 
ratio, or lumi nosity. Bas e d on this and high luminosities in both the i^-band and [OIII] 
emission line, iDong et al.l (120051 ) argued that objects with steep Balmer decrements in their 
sample are intrinsically luminous objects reddened by dust. On the other hand, steepening 
of Ba lmer decrement when c ontinuum dimmi n g has been repor ted for individual objects 
(e.g., IWamsteker et al.l Il990l ; iTran et al.l Il992l ). iGoodrichI (119891 ) argued that the variable 



Balmer decrement in some objects of their sample can be attributed to dusty clouds moving 
into/out of the line of sight. 

The variable Balmer decrement in Mrk 1393 can be most naturally interpreted as vari- 
able extinction. First, X-ray absorption with a column density around 2 4x10^^ cm~^ by 
cold material, as expected, has been detected in this source when it showed a steep Balmer 
decrement of Ha/H(3 ~ 6.0. Second, the observed broad Ha flux during the BAO obser- 
vation in 2005-2006 with a steep Balmer decrement is very similar to that in 1984 when 
the Balmer decrement was normal, indicating that the steep Balmer decrement detected by 
BAO is not caused by continuum luminosity variations as observed during the lower state of 
NGC 5548. Also consistent with continuum extinction is the fact that Mrk 1393 displayed 
a blue continuum in 1984 when the Balmer decrement is normal, while all of the reddened 
spectra show a stellar dominated continuum. We note that the stellar dominated continuum 
in 1993 is likely due to the intrinsic weakness of the continuum, since the Ha broad line flux 
is about 4 times lower than in 1984. 



Now that it has been established that the steep Balmer decrement in Mrk 1393 is caused 
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by dust reddening, we can use the Balmer decrement to correct for the dust extinction. The 
extinction corrected broad Ha flux was the lowest in 1993 (at the KPNO observation), and 
the highest in the end of 2005 and early 2006 (during our BAO observations), giving the 
largest observed variation of a factor of 7.5. The intrinsic Ha flux was in between these values 
at the epoch of the 1984 observation iMorris &: WardI (119881 ) and the 2001 SDSS observation 
(a factor of 1.9 and 2.1 higher than the lowest state, respectively). These broad Ha line 
flux variations do not correlate with the Balmer decrements. X-ray flux during the ROSAT 
observation in 1991 is a factor of six lower than in the end 2005 and early 2006, perhaps 
indicating the object was in a low state in the early 1990s. We suggest that Mrk 1393 is 
a good candidate for long-term monitoring. Because the Balmer decrement varies, as what 
was seen from 1993 though 2001 to 2006, the object could be an ideal laboratory for detailed 
studies of the physical and geometric properties of obscuring matter. 



4.2. Dust-to-Gas Ratio of the obscuring gas 



Since the BAO spectra were taken only one and half month after the XMM-Newton 
observation and nearly simultaneous with the Swift observation, it is reasonable to assume 
little or no variation for the BLR reddening on such short time scales. The quasi-simultaneity 
of the X-ray and optical observations enables us to derive the dust-to-gas ratio for the 
obscurin g gas in Mrk 1393 . Assuming an intrinsic Ha/H/3=3.0, the average value for blue 
quasars JPong et al.ll2008h . Mrk 1393 was reddened by E{B - l/)=0.59 during the BAO 
observations. The column densities derived from the X-ray absorption from the XMM- 
Newton and Swift data, are around 2.0-3.6 xlO^^cm"^ for different spectral models assuming 
solar abundances. We find a dust to gas ratio of E{B — V)/Nh = 1.6 — 3.0 x 10"^^ mag cm^ 
if the absorber is cold. In the two phase ionized absorber model for X-ray spectrum, the 
column density of the low ionization material was found to be 5.4xl0^^cm~^, which gives 
a dust to gas ratio of 1.1x10"^^ mag cm~^. These values are within a factor of two of the 



average Galactic value 1.7x10"^^ mag cm~^f lBohlin et al.lll978l ). Furthermore, the strong 
variability in red dening in Mrk 1393, as well as in several other intermediate Seyfert galaxies 
( lGoodrichlll989l ) . warns against the use of non-simultaneous optical and X-ray data to derive 
the dust-to-gas ratio, at least for individuals. 

As a comparison, the observed dust-to-gas ratios in other AGN show a large diversity. 
A class of (low- and high-luminosity) AGN has emerged which show relatively high optical- 
UV extinction, but lack a corresponding amount of (cold) X-ray absorption (e.g., Brandt 
et al. 1997, Komossa & Fink 1997b, Leighly et al. 1997, Komossa & Bade 1998). Three 
mechanisms have been proposed to account for these observations: (i) a higher-than-Galactic 
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dust/gas ratio, (ii) a different grain-size distribution with a dominance of small grains, and 
(iii) dust mixed with absorbing material which is ionized rather than neutral (Komossa & 
Fink 1997b). In at least some 10 AGN the third solution is favored (see Komossa 1999 for a 
review). On the other hand, Maiolino et al. (2001) presented a sample of mostly luminous 
Seyfert galaxies which show indications for dust /gas ratios much lower than the Galactic 
value (by a factor 3-100). It is possible that the observed diversity of the dust to gas ratios 
in AGN can be reconciled by a combination of different obscuring components. 



4.3. On the Nature of the Obscuring Gas 

The obscuring material seen in Mrk 1393 must be very close to th e AGN. By a coni- 



parison of the Balmer decrements of the broad and narrow components, iDong et al.l (120051 ) 
pointed out that the reddening of the NLR is much smaller than that of the BLR in most 
of the partially obscured quasars in their sample. The authors then suggested that the 
obscuring material lies in between the BLR and the NLR. The observed variations of the 
reddening set constraints to the size to an even smaller region in Mrk 1393. The object was 
not seriously reddened in 1993 (the KPNO observation) but was heavily reddened in 2001 
(the SDSS observation) and 2005-2006 (the BAO observations). We interpret this as the 
obscuring matter moving into (~1993-2001) and out of (~2001-2006) the line-of-sight to the 
nucleus (both of the continuum source and the BLR), as dust condensation likely takes a 
much longer time in the vicinity of the galactic nucleus. The size of the BLR, estimated from 
the Ha line luminosity (reddening correct ed) and the emp i rical Rrt,r versus Ha luminosity 



relationship, is R ~15 light days, following I Wang fc Zhang] (120031 ). Using the revised relation 



given by Kaspi et al. (2006) results R ~19 light days. The obscuring material must cover 
most of the BLR, and thus the transverse velocity of the obscuring material must be not less 
than 1000 km s~^. With this transverse velocity, the absorber is unlikely much farther than 
a few parsecs away from the nucleus. As such, it may be part of the putative torus or even 
within it. If this is the case, our result suggests that a dust-to-g as ratio simil ar to the Galac- 



tic value can be found for gas located very close to the AGN. iMattI (120001 ) proposed that 
Compton thin Seyfert 2/1.8/1.9 galaxies are absorbed by dust lanes in their host galaxies. 
However, this is not the case for Mrk 1393. 

Certain constraints on the gas properties can be inferred from the data. Assuming that 
the reddening variability is caused by a single cloud event and the cloud is quasi-spherical, 
we can set an upper-limit on the particle density of the cloud of the order of 10^ cm~^, 
by using the X-ray absorption column density and a cloud size no less than the BLR size. 
Alternatively, if the cloud is pancake-like under the radiation pressure, the density of the gas 
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can be higher. The cloud is unhkely self-gravitational bound because the tidal force of the 
central black hole exceeds the self-gravity of the cloud by 1-2 orders of magnitudes for the 
above estimated density, if the distance to the black hole is less than 10 pc. The disrupted 
cloud will have a filament structure and the low column density gas will be blown away due 
the strong AGN radiation pressure. Thus the material must be either transient or bounded 
by some external pressure against tidal disruption. 

Instead of appealing to the unknown external pressure, the transient scenario is much 
more natural. One natural possibility is that the absorbing gas is debris disrupted from 
the dusty torus. The black hole mass of Mrk 13 93 is estimated to be 5.5 xlO ^M© from 
the Mhh — (J relation]^ (e.g. Tremaine et al. 2002 ). and 2x10® from the above BLR size 



and the emission line width {ana = 2940 km s~^). The bolometric luminosity estimated 
from the observed hard X-ray luminosity (9x10^^ ergs s~^) is about 0.10—0.37 of the Ed- 
din gton luminosity, assu ming a bolometric correction to the 2-10 keV X-ray luminosity of 



30 (IMarconi et al.ll2004l ). Fabian et al. (2008) showed that the radiation force exceeds the 
gravity of a black hole for a dusty gas below the column density Nh < 5 x lO^^A cm~^, when 
A is the Eddington ratio. With A = 0.10 — 0.37, the measured X-ray column density in Mrk 
1393 is much lower than this critical value, and thus the gas must be accelerated outward 
by the radiation force if the black hole is the only source of gravitational force. The masses 
of gas and stars within a few parsecs do not likely exceed that of the black hole by a large 
factor, so they are insufficient to prevent the outflow. Thus, the obscuring gas we are seeing 
is either formed in the inner nuclear region, or bound and shielded by thick material before 
it is exposed to the nuclear radiation. However, in the presence of a strong radiation fleld, 
it is unrealistic to form dust in relatively thin material within the inner parsecs. A plausible 
origin of the gas is thus debris disrupted from the dusty torus. 

To summarize, our XMM-Newton and Swift observations detected signiflcant X-ray ab- 
sorption from cold gas with a column density of about 2 — 4 x 10^^ cm~^ in the partially 
obscured luminous Seyfert galaxy Mrk 1393. Our semi-simultaneous BAO optical spectro- 
scopic observations gave a steep Balmer decrement of Ha/HjS ~ 6.0. A combination of 
the X-ray and optical observations suggests that the steep Balmer decrement is caused by 
dust reddening, and the dust-to-gas ratio of the obscuring material is similar to the Galactic 
value. Incorporating the archived and literature data, we flnd variations in Balmer decre- 
ment on a time-scale of several years. Such variations require that the obscuring material 
locates within a few parsecs from the central nucleus. We suggest that the obscuring material 
may be debris evaporated/disrupted from the dusty torus. The source also displayed large 



6 



We measured a stellar velocity dispersion of 159±5 km s ^ in the SDSS spectrum following the method 



of ILu et al.l (|2006f ). which is consistent with 163±15 km s ^ as given by Greene & Ho (2006). 
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amplitude X-ray variations of a factor 6 over 15 years. 
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Fig. 1. — Left panel: The best fit to the XMM-Newton spectrum in the 0.3-10 keV band 
with a model consisting of a broken power-law (short-dashed) absorbed by cold gas and 
thermal plasma emission (dotted) with the PN data in green and MOS data in red and 
black. The individual components of the best-fit model are also shown. Right panels: 
residuals for different models (from top to down): a) powerlaw model extrapolated from 2- 
10 keV; (b) partial covering absorption model; (c) single ionized absorption model (absori); 
(d) absorbed powerlaw plus meka model; (e) power-law-|-black body emission absorbed by 
cold gas plusmeA;a model; (f) broken power-law absorbed by cold gas plus meka model 
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Fig. 2.— The contours at 70%, 90% and 95% (Ax^ =2.3, 4.6 and 9.2) confidence levels on 
the plane of the intrinsic absorption column density A'^^ versus photon index for model-c 
(red) and model-d (black) in Fig 1. Model-e gave a contour very similar to model-d, thus 
not show here for clarity. 
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Fig. 3. — The background-subtracted X-ray light curve during Swift XRT observation. The 
curve is rebined to have 100 counts per bin. Short time scale variability is apparent. 
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Fig. 4. — The optical spectra of Mrk 1393 taken at three epochs are normahzed to [OIII] A5007 
and shown in the upper-left panel (see text for details). The KPNO spectrum is offset by 
20 X 10~^^ erg cm'"^ for clarity. The lower-left panel is the same plot expanded in 
the vertical axis for a clear view of the continuum and broad lines. Our best fit continuum 
model is also over-plotted in the panel. Emission line spectrum (after subtraction of the 
best fit continuum model) is shown in the right panels (upper and middle panels for the 
H/3+[0III] regime, and lower panel for the [OI]+Ha+[NII] + [SII] regime). It is clear that 
the continuum is dominated by starlight in all of the KPNO, SDSS, and BAO spectra. Also 
note the dramatic variability of the broad Ha and H/5 lines, while all of the narrow lines 
remain constant in the three spectra. 
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Fig. 5. — The overall variation of the broad Ha flux ( "0" ) and the reddening ( "□" ) estimated 
from the Balmer decrement. 
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Table 1. Results of acceptable X-ray spectral fits 



model 


Nh 




r 


F2-10keV kT''''/r2 




-'^ 0.5-2. OkeV 




(1) 


(2) 




(3) 


(4) (5) 


(6) 


(7) 


(8) 










XMM-Newton 1-10 keV 






1 


9 9+0.56 
^•^^-0.60 


1 


•^"-"-0.06 


8.17 






431/478 










XMM-Newton 0.3-10 keV 






2 


q nO+0-46 
O.UZ_Q 47 


1 


cq+O.OS 
•^'-'0.05 


8.12 


Q if;f;+0.004 
U.iUU_Q 004 


5.69 


667/650 


3 


9 07+0.04 
^•'5' -0.03 


1 


/I 8+0-05 


8.21 0.054+°:°°^ 


0.167^^°^ 


3.59 


643/649 


4 


9 9 f- +0.03 
■^•■^^-0.04 


1 


/1S+0-04 


8.10 5.69^?;^^ 
Swift 1-10 keV 


0.1681^?^ 


3.09 


634/650 


1 


0+1.1 
O.O_2.0 


1 


c-o+0.21 
•^0-0.19 


7.5 

Swift 0.5-10 keV 






34/51 


4 


3.611^ 


1 


Rn+0.2A 
•""-'-0.21 


7.5 4.51^1 


1 c^fi+0 041 
U.iOO_o 078 


3.18 


43/62 



^model: Galactic absorption has been applied to all models. The absorption column 
density listed in this table is referred to the intrinsic absorption at the redshift of Mkn 1393. 

1 - a power-law absorbed by cold gas in the source rest frame. 

2 - a power-law absorbed by cold gas plus thermal plasma emission {meka in Xspec) 3 - 
a power-law and a black body absorbed by cold gas plus thermal plasma emission [wabs x 
{zwabs{powerlaw + zbb) + meka) in XSPEC\. 

4 - a broken power-law absorbed by cold gas plus thermal plasma emission [wabs{zwabs x 
bknpow + meka)] 

^{2) absorption column densities in units of 10^^ cm~^; (4) modeled X-ray flux in units of 
10~^^ ergs cm~^ s~^; (5) kT of black-body component in units of keV; (6) kT of meka compo- 
nent in units of keV; (7) flux of meA;a component in 0.5-2.0 keV in units of 10~^^ ergs cm~^ s~^. 



